Zigzag spin chains with antiferromagnetic-ferromagnetic interactions: Transfer-matrix 

renormalization group study 
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Properties of the zigzag spin chains with various nearest-neighbor and next-nearest-neighbor in- 
teractions are studied by making use of the transfer-matrix renormalization group method. Ther- 
modynamic quantities of the systems (temperature dependence of the susceptibility and the specific 
heat), as well as the field dependence of the magnetization are analyzed numerically with a high 
accuracy in the thermodynamic limit. The results have been compared with the recent experimental 
data on Rb2Cu2Mo30i2- 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Mg 
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I. INTRODUCTION 

Spin models with charge degrees of freedom frozen, as 
an effective low-energy description for insulating systems, 
have been proved to be a fruitful resource of fundamen- 
tal concepts and principles in condensed matter physics 
and other related fields. Besides theoretical explorations, 
advanced probe tools and measurements, combined with 
multifarious natural and artificial materials, such as ce- 
ramic or organic compounds, and optical lattices provide 
remarkable practical platforms for theoretical investiga- 
tions. The simplest one-dimensional (ID) quantum spin 
model is the nearest-neighbor (NN) Heisenberg model: 
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where Si is the quantum spin operator. From Bethe's 
seminal papesi and successive works, we know that the 
ground state of a spin-i antiferromagnetic (AF) (J > 0) 
chain is a singlet (S to t = 0) and has quasi-long-range or- 
der with algebraically decaying spin correlations. The 
gapless spectrum contains no single-particle excitations 
and is instead a continuum of states. The elementary 
excitations are called spinons, which carry spin 1/2 and 
appear only in pairs in all physical states with integer 
total spin. This picture is qualitatively different from 
the prediction of a spin-wave theory, which is usually 
effective in higher-dimensional systems. The gapless be- 
havior is special for half-integer Heisenberg spin chains, 
while for integer ones, Haldane conjectured that there 
exists a finite gap^ Experimental, numerical, and theo- 
retical studies have confirmed this conjecture for S = 1 
and some other higher spin values^ On the other hand, 
half-integer spin chains can be driven to a gapped phase 
by frustrations. Throughout this paper, we confine our 
discussions to the spin- 1/2 case. 

A straightforward generalization of the NN Heisenberg 
model is to include the next-nearest-neighbor (NNN) in- 
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FIG. 1: Heisenberg zigzag spin chains. 
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where J\ or J 2 is NN or NNN exchange interaction 
constant, and periodic boundary conditions are implied. 
This is usually called the J1-J2 model, which can be also 
considered as a model for a zigzag spin chain, as shown 
in Fig. HJ 

The J1-J2 model has been investigated theoretically 
over the decades. With J 2 > 0, it is a frustrated (com- 
peting) system, irrespective of the sign of J\. When 
J\ > and J 2 > (AF-AF), the ground state is a 
spin liquid. The increase of the ratio of the coupling 
constants a(= J2/J1) induces an infinite-order phase 
transition from a gapless state to a gapful dimerized 
stateA^ The critical point a c is numerically estimated to 
be sw 0.241^ When a is further increased to the so-called 
Majumdar-Ghosh (MG) point at 1/2, the ground state is 
the products of singlet pairs formed by nearest neighbor- 
ing spinsi^ It is two-fold degenerate since the Z2 sym- 
metry of translations by one site is spontaneously broken. 
When Ji < and J 2 > (F-AF) with -1/4 < a < 0, 
the ground state is fully ferromagnetic (FM), and be- 
comes an (S = 0) incommensurate state for a < — 1/4. 9 - 
It is suggested that in this incommensurate state the gap 
is strongly suppressed^ For a = —1/4, the exact ground 
state can be shownii to have a (N + 2)-fold degeneracy, 
comprised by S = and S = N/2 states (with N the lat- 
tice size). When J\ > and J 2 < (AF-F), the system 
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is believed to be in a gaplcss antifcrromagnctic phase for 
any permissible values of J\ and J 2 . 

There are many papers contributing to the AF-AF case 
and the related extended models. The other two cases (F- 
AF and AF-F), although having caught relatively less at- 
tention, also show very interesting phenomena, which we 
will mainly deal with in this paper. Aside from the gen- 
eral aspect of theoretical interest, especially for under- 
standing the roles played by frustration and incommen- 
surability, an additional motivation to study this system 
lies in the fact that physically, it is believed that a large 
class of copper oxides can be essentially described by the 
J1-J2 model. 

Copper oxides are excellent model systems for low 
dimensional spin-^ quantum magnets, where magnetic 
Cu 2+ ions carry 1/2 spins. The basic building blocks are 
CuC-4 plaquettes and there are three ways of linking these 
fundamental units. One is adjacent squares sharing their 
corners, as shown in (a) of Fig. A typical example is 
Sr2Cu03. These corner-sharing chains can expand in the 
plane to form a Cu02 sheet, constituting the basic struc- 
ture in cuprates. The dominant interaction in a corner- 
sharing chain is the NN superexchange. Linear Cu-O-Cu 
bonds along the spin chains give rise to a large antifer- 
romagnctic NN exchange coupling. As in Sr 2 Cu03, the 
NN coupling constant is estimated to be 2100 ± 200 
The second kind is edge-sharing, as shown in Fig. 0(b). 
Because of the nearly 90° Cu-O-Cu bond in the edge- 
sharing squares, an 02p a orbital hybridizing with a 3d 
orbital of Cu ion is almost orthogonal to that of the next 
Cu ion. The NN interaction J\ can vary from antiferro- 
magnetic to ferromagnetic, as the angle 9 of the Cu-O-Cu 
bond approaches 90° from a larger valued This nearly 
orthogonality makes the NN coupling in the edge-sharing 
case more than an order of magnitude smaller than the 
corner-sharing. The sign and the absolute value of the 
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FIG. 2: Various ways of linking CUO4 plaquettes. (a) 
The corner- sharing chain, as in Sr2Cu03; (b) the edge- 
sharing chain, as in CuGeOa, Li2Cu02, and Rb2Cu2Mo30i2; 
(c) Cu02 plane expanded by corner-sharing chains, as in 
cuprates; (d) the zigzag chain in SrCu02- 



NN interaction depend sensitively on the bond angle 6 
and the distance between copper and oxygen ions. We 
refer readers to Table I of Ref. ^| for details. The third 
configuration is the combination of corner-sharing and 
edge-sharing in one spin chain simultaneously, such as in 
SrCu0 2 , as shown in Fig. 0(d). In SrCu0 2 , the NNN 
coupling, which is the superexchange interaction through 
the linear Cu-O-Cu, is almost ten times greater than the 
NN one through sharing the cdges>i£ 

In corner-sharing chains, the NNN interaction can of- 
ten be neglected safely due to its small magnitude rel- 
ative to the NN interaction. For the edge-sharing case, 
the situation is quite different. The NNN interaction J 2 
through the Cu-O-O-Cu path is generally antiferromag- 
netic, usually with a magnitude of a few tens Kelvin. 
Despite its smallncss, J 2 has a pronounced effect on the 
physical properties of these systems since the NN cou- 
pling is also small. Therefore, edge-sharing copper ox- 
ides provide abundant experimental materials for study- 
ing the zigzag spin chain model, which can cover a large 
region of the parameter space. 

On the other hand, at low temperatures, besides in- 
trachain couplings J\ and J 2 , other forms of interac- 
tions often become relevant, driving a system to vari- 
ous phases with the decrease of temperature. For exam- 
ple, spin-phonon interactions can induce a spin-Peierls 
instability, as in CuGeOa i& If interchain interactions are 
strong enough, an antiferromagnetically long-ranged or- 
der will appear below the Neel temperature T/v, as in 
Ca 2 Y 2 Cu 5 10 £L La 6 Ca 8 Cu 24 4 i and Li 2 Cu0 2 Mi An- 
other interesting case is LiCu 2 2 j2£ which undergoes a 
transition to a magnetic helix state at low temperatures. 
By a comparative study on the similar oxide NaCu 2 2 , 
it has been shown that interchain interactions over a few 
chains should be incorporated to explain the experimen- 
tal results^ These facts reveal the complexity of the in- 
teractions underlying edge-sharing copper oxides at low 
temperatures. 

Since the effect of NNN interaction is important for 
edge-sharing copper oxide chains, the study on the J\- 
J 2 model not only has its own theoretical meaning, but 
also obtains a connection with practical materials. In 
this paper, in the thermodynamic limit and extending to 
the low-tempcrature region, we use the transfer-matrix 
renormalization group (TMRG) methooj22*2ii24 to study 
the thermodynamic properties of the J1-J2 model with 
antifcrromagnetic- ferromagnetic interactions. 

The TMRG is a powerful numerical tool for studying 
the thermodynamic properties of ID quantum systems. 
It starts by expressing the partition function as a trace on 
the product of the transfer matrix Tm using the Trotter- 
Suzuki decomposition. 



Z = Tre-P H = lim TrTj /2 , 



(3) 



where M is the Trotter number, and TV is the total cell 
number in the lattice, r =_0JM. The definition of Tm 
can be found from Refs. S22I2I2I For the J x -J 2 model 
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considered here, each cell consists of two adjoining spins. 
In the thermodynamic limit TV — ► oo, the free energy 
/. internal energy u, and uniform magnetization m z can 
be expressed by the maximum eigenvalue Amax and the 
corresponding left (tp L | and right | ip R ) eigenvectors of 
the transfer matrix Tm- 
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where the definition of the transfer matrices Tjj and Tm , 
which are similar to Tm, can also be found in Refs. I22I 
I23TI241 The specific heat and magnetic susceptibility can 
then be calculated by numerical derivatives of u and m z , 
respectively, 
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In our numerical simulation, r = 0.05, the error caused 
by the Trotter-Suzuki decomposition is less than 10~ 3 . 
During the TMRG iterations, 60 — 80 states are retained 
and the truncation error is less than 10~ 4 down to ksT ~ 
0.01J. 

The outline of this paper is as follows: Section [n] 
and Sec. Mil are devoted to discussions on F-AF and 
AF-F cases, respectively. Section ITVl shows our numer- 
ical results compared with the newly experiments on 
Rb2Cu2Mo30i2^ We conclude with a brief summary 
in Sec. 



II. THE F-AF CASE 

As mentioned in the previous section, at the critical 
point a c = —1/4, two distinct configurations with the 
energy E g = — 3TV|Ji|/16 are the ground states, where 
A is the lattice sizeiii One is fully ferromagnetic with 
Stat = A/2, the other is a singlet state with S tot = 0. 
The state vector for the latter can be expressed as 



® = J2[iJ] [M [m,n] 



where the summation is made for any combination of 
spin sites under the condition that i < j, k < I, m < n, 
■ ■ •, and denotes the singlet pair. This is also called 
a uniformly distributed resonant- valence-bond (UDRVB) 
state. 

In the region 0>a>— 1/4, the ground state lies in the 
subspace Stot = A/2 with the degeneracy A + 1. The 
ground state energy E g = —N\J\\ (1 + a) /4. For a < 
— 1/4, the ground state lies in the subspace Stot — Stot = 



and the lattice translational symmetry is thought to be 
broken. For the critical point a c = —1/4, besides the 
ferromagnetic configuration, the UDRVB state restores 
the lattice translational symmetry. 

When a < —1/4, whether the system is gapped or 
gapless is an interesting and controversial issue. When 
Ji \Ji\, it is appropriate to regard the model as two 
antifcrromagnetic spin chains coupled by a weak zigzag 
interchain interaction J\ . This coupling can be expressed 
by the current-current interaction^^ in terms of the 
Wess-Zumino-Witten fields (see, e.g., Ref . l26h . If J\ > 0, 
by renormalization group (RG) analysis, this interac- 
tion is marginally relevant and produces an exponentially 
small gap A cx exp (— const J2/J1), which leads to a spon- 
taneously dimerized ground state£2ii While at the ferro- 
magnetic side, i.e. J\ < 0, the model was believed to be 
gapless because the current-current interaction renormal- 
izes logarithmically to zero^ It was conjectured that the 
ferromagnetic model is critical with different velocities for 
the spin-singlet and spin-triplet excitations^ Afterward, 
Nersesyan et alSL found that in addition to the current- 
current interaction, a "twist" term associated with the 
staggered component of the spin operators arises in the 
zigzag chains. Due to this parity-breaking term, the crit- 
ical point Ji = is unstable both in the ferromagnetic 
(Ji < 0) and antifcrromagnetic (Ji > 0) regionsi 27 i 28 

The phase diagram for the F-AF model in an external 
magnetic field was discussed by ChubukoviS In addition 
to the ferromagnetic phase, two different biaxial and uni- 
axial spin nematic phases are mapped out. In these ne- 
maticlike phases there is an extra symmetry breaking of 
reflections about a bond or about a site. In the absence of 
external field, with the decrease of a starting from —1/4, 
the system develops from the chiral biaxial spin nematic 
phase to the dimerized uniaxial spin nematic phase at 
a ~ —0.385. The Licb-Schultz-Mattis theorem^ states 
that a half-integer spin chain with essentially any rea- 
sonably local Hamiltonian respecting translational and 
rotational symmetries either has a zero gap (i.e. "mass") 
or else has degenerate ground states, corresponding to 
a spontaneously broken parity. If the phase diagram is 
correct, we could expect an energy gap at a < —1/4. 

On the other hand, numerical analysis has shown a 
complicated size dependence of the ground-state energy 
and correlation function in the region a < — 1/4 £2i This 
phenomenon, combined with the fact of slow convergence 
and no detection of energy gap at the resolution of the 
numerical simulations indicates an unusually long corre- 
lation length in the F-AF chain. 

By taking into account the twist term, the RG 
analysis— found that although the unstable RG flow 
around the critical point J\ = produces an energy gap 
in the ferromagnetic coupling as well as the antiferromag- 
netic one, in the ferromagnetic side, due to the existence 
of a marginally relevant fixed line, the gap is strongly re- 
duced. In an extended region for the J\ of order one, the 
correlation length can be extremely large and the gap, if 
exists, is so strongly suppressed that numerical methods 
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can not detect it. 

Another closely related topic worth mentioning is the 
incommensurability in the zigzag spin chains. For sim- 
ple antiferromagnetic Hciscnberg spin-1/2 chains and 
ladders, the well known mechanisms for generating in- 
commensurabilities are via external magnetic fields or 
Dzyaloshinskii-Moria interaction. While in the zigzag 
spin chain, it has been found that the frustrated interac- 
tion J 2 > can also produce incommensurability^i 

In classical picture, by regarding spin operators Si as 
classical vectors, the energy per site can be expressed as 

E (9) = Ji cos 6 + J 2 cos 26> (9) 

for Hamiltonian @. The pitch angle is given by 

cos (9 = -Ji/U 2 (10) 

for |a| = I J2/J1I > 1/4. As Ji > 4J 2 > 0, 9 = tt; -Ji > 
4J 2 > 0, 9 = 0. At the special point J\ = 0, 9 = ir/2, and 
the deviations from the normal values (tt for the AF-AF 
case and for the F-AF) begin to occur at |a| = 1/4. 
In quantum level, the characteristic momentum Q of the 
spin-spin correlation function which maximizes the static 
structure factor is either tt or 0, corresponding to the 
antiferromagnetic or ferromagnetic spin chain with J 2 
0. respectively. With the increase of J 2 , Q departs from 
its usual value. Numerical simulations^^ display that 
for the AF-AF case, Q deviates from tt with the increase 
of J 2 after crossing the MG point (J 2 = Ji/2) where the 
dimerization has already taken place. On the other side, 
the departure of Q from sets in at the critical point 
J 2 = — Ji/4 for the F-AF case. This asymmetry may 
be understood on account of the effects of the current- 
current interaction and the twist term. 

It was shown that in the presence of exchange 
anisotropic, the twist term induces incommensurabilities 
in the spin correlations^ We have reason to expect this 
to hold true even in the SU(2) symmetric (i.e. isotropic) 
casei 32 i 33 In the AF-AF isotropic case, both the twist 
term and the current-current interaction diverge simul- 
taneously to reach the strong coupling phase as the RG 
flows to J 2 3> J\ ~ The pure current-current interaction 
induces massive spinons which can be regarded as quan- 
tum dimerization kinks, driving the zigzag spin chain into 
the dimerized phase. The twist term appearing in the 
zigzag case is merely to shift the minimum of the two- 
spinon continua to incommensurate wave numbers, which 
does not alter the massive spinon picture qualitatively. 32 
If these two terms become relevant at different points, 
we may observe the emergence of dimerization and in- 
commensurability one after the other, as in the AF-AF 
case. For the F-AF case, the current-current interaction 
is marginally irrelevant since it renormalizes logarithmi- 
cally to zero^ The twist term becomes dominative. We 
can observe the incommensurability only, and hardly de- 
tect the gap following with the dimerization. But there 
is not a simple way to separate the effects of the current- 
current and twist interactions in the isotropic zigzag spin 
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FIG. 3: Temperature dependence of the uniform susceptibility 
at various a for the F-AF case. 



chains as in the anisotropic ones. There is still much the- 
oretical work to be accomplished. 

Since the TMRG gives results for observable quantities 
in the thermodynamic limit, we are able to exploit it to 
study the bulk properties of the system without worrying 
the finite size effect. The complicated size dependence of 
the ground-state energy and correlation function found 
in the region a < —1/4 (RefsA^i) may be avoided. 

Figures|31and0]show the TMRG results on the temper- 
ature dependence of the susceptibility Xi specific heat C, 
and heat coefficient C /T at various a for the F-AF case. 
At a = 0, the model reduces to the case of ferromagnetic 
spin chain. The curves of a = describe the properties 
of x an d C for a ferromagnetic spin chain, which behave 
in the low temperature limit: 

X ~ T-\ 

c ~ VT. 

At a > —1/4, the temperature dependence of x always 
diverges, indicating that the system lies in a ferromag- 
netic state. At the critical point a = —1/4, a phase tran- 
sition from the ferromagnetic state to the singlet state 
takes place. In the region a < —1/4, remarkable sup- 
pressions of the susceptibility can be observed, and with 
the decrease of a, the peaks of x move to higher temper- 
atures with its heights decreased rapidly. 

For the temperature dependence of the specific heat, 
the most remarkable feature is the development of a 
double-peak structure. With the decrease of a, at inter- 
mediate temperatures, a broad maximum is maintained 
and its height lowered. A relatively sharp peak at low 
temperatures, induced by the NNN AF interaction J 2 , 
appears and develops, and its position moves to higher 
temperatures. We expect that this peak would approach 
the maximum of C(T) for a pure AF Heisenbcrg chain^4 
The drastic change on the shape of the curves between 
a = —0.125 and a = —0.28 implies a phase transition. 
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FIG. 4: The specific heat C and heat coefficient C/T at var- 
ious a for the F-AF case. 



The similar double-peak structure for the Ji-J? model 
in the F-AF case has also been found in Ref. |35j By 
applying a method of hierarchy of algebras, the authors 
calculated thermodynamic quantities of the linear ring 
of size 16 described by the same model (|2J). Our results 
arc qualitatively in agreement with theirs. Recently, this 
double-peak structure in the specific heat for the mod- 
erate value a = —1/3 has also been confirmed by ex- 
act diagonalization methods^ Furthermore, for various 
a we have considered, C/T as a function of T decreases 
monotonously down to T ~ 0.03. This fact reflects a high 
density of low-excitation states in this region, and con- 
firms the results obtained from the density-matrix renor- 
malization group calculations^ 

From the numerical results above, we find that the ex- 
istence of gap near the critical point a c is still a question. 
The behaviors of susceptibility and specific heat indicate 
that if a gap exists near a < —1/4, it should be very 
small. The properties of the excitations near and far 
from a c need to be investigated further. 



III. THE AF-F CASE 

In the AF-F case with Ji > and J2 < 0, the behavior 
of the zigzag spin chain is believed to be antiferromag- 
netic with no gap for any permissible values of J\ and 
J 2 ^i The arguments arc based on a simple spin- wave 
analysis. The spectrum of spin- wave excitations is given 
by ' 



e k = 2SJ\l 



where S = 1/2, a = J2/J1 < 0, X k 
7fc = cos k. efc is linear as k — > 0: 



11] 



-2a sin k + 1. 



vk 



(12) 



with v — 2S\/—Aa + 1. Although the spin- wave theory 
is not quite correct for antifcrromagnctic spin chains, it is 
still insightful in shedding light on the qualitative picture 
of the excitations around the characteristic momentum 
k = 0, 7r. We see that for the AF-F case, there also exist 
gapless excitations at k = 0, n. 

According to Lieb's discussion^ Hamiltonian @ with 
Ji > 0, J2 < describes spin systems at a bipartite lat- 
tice. The absolute ground state of the AF-F spin chain 
should lie in the S = sector. On the other hand, the 
relative ground state in each subspace V(M) of the to- 
tal magnetic quantum number M is unique. Thus we 
deduce that the absolute ground state in the AF-F case 
is nondegenerate. The Lieb-Schultz-Mattis theorem as- 
serts that for a half-integer spin chain with reasonably 
local Hamiltonian respecting translational and rotational 
symmetries, either the parity is spontaneously broken in 
the ground state or else there are gapless excitations of 
odd parity, under the condition that the ground state is 
rotational invariant^ Combining these two statements, 
we expect the system should be gapless. 

It is not trivial to make the above statement completely 
rigorous, considering the relative ground state in the sub- 
spaces V(M) may not be unique in the thermodynamic 
limit N — > 00. In order to confirm the conclusion, wc 
calculate the temperature dependence of the susceptibil- 
ity and specific heat for various a by using the TMRG 
method. 

Figures|5]and||)]show the TMRG results on the temper- 
ature dependence of the susceptibility x, specific heat C, 
and heat coefficient C/T at various a for the AF-F case. 
With the decrease of a, the behaviors of x an d C do not 
change qualitatively, x — * const (nonzero) and C ~ T as 
T — > 0. The Bonner-Fisher peak in x(T) (Refi2£) moves 
to higher T with the decrease of a. Phase transition 
is not observed. The slopes of C, or the intersections 
of C/T at zero temperature (shown in Fig. 0) decrease 
with the decrease of a. This indicates that the density of 
low-excitation states becomes smaller and the spin fluc- 
tuations are suppressed by the NNN FM interaction. 
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J 1 > 0, J 2 = aJ r a < J 1 > 0, J 2 = aJ v a < 




FIG. 5: Temperature dependence of the uniform susceptibility 
at various a for the AF-F case. 



IV. COMPARISON WITH THE 
EXPERIMENTAL RESULTS 

In the previous two sections, we have discussed ther- 
modynamic properties of the J\-J% model for the F-AF 
and AF-F cases and our numerical results cover their 
whole phases. In experiments, in order to understand 
the physics of edge-sharing copper chains, it is impor- 
tant to determine both sign and magnitude of NN and 
NNN interaction coefficients. One way for this pur- 
pose is to measure thermodynamic quantities, such as 
specific heat, susceptibility and magnetization, together 
with numerical fitting on these data. Such a method 
has been widely used in the studies on copper chains 
and its efficiency has been proved, such as in SrCu02r*" 
NaCu 2 3? 2i Ca2Y2Cu.5O1.5jii and La 6 Ca 8 Cu 24 04i^ 

Recently, the magnetic susceptibility and mag- 
netization of the edge-sharing copper oxide 
Rb 2 Cu2Mo 3 Oi2 (Refi2&) have been measured.^ The 
most interesting finding is that no magnetic phase 
transition was observed down to 2 K. Therefore, the 
compound is suitable for studying the properties of 
the ground state of the J\-Ji model. It was proposed 
that at first approximation, it could be described by 
a F-AF J1-J2 models We have known that for the 
F-FA chain, the unusually long correlation length, which 
exist even at high-frustrated region, can lead to promi- 
nent finite-size effect for the usual cluster simulations. 
Therefore, as a primary step, reliable numerical results 
on thermodynamic quantities free of finite-size effect 
on the pure J1-J2 model may be essential in making 
a quantitative comparison between experimental and 
theoretical studies. In this section, we use the TMRG 
method to simulate the experimental results based on 
model ©. 

The strategy is as follows: For a fixed a, we first use the 
TMRG to calculate the temperature dependence of the 
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FIG. 6: The specific heat C and heat coefficient C/T at var- 
ious a for the AF-F case. 



susceptibility in the units of J±. Then, according to the 
position of the peak in the numerical result and the actual 
value obtained in the experiment, we can determined the 
value Ji uniquely, on the condition that the calculated 
peak's position coincides with the exp eriment. The g 
value is taken to be 2.03 as in Ref. ITiL The following is 
what we observed. 

For the AF-AF case (Fig. 0), the height of the peak 
decreases with the increase of a. The possibility for this 
type of interaction can be excluded. 

For the F-AF case (Fig. [SJ), we see that with a ap- 
proaching the critical point —1/4, to keep the peak po- 
sition unchanged, the magnitude of J± acquires a value 
of thousands of Kelvin. Furthermore, the deviation from 
the experiment data at high temperatures becomes more 
conspicuous. While taking the parameters J\ = —138 K 
and a = —0.37 as in Ref.[lJ, the temperature for the peak 
is less than T max = 14.3 K found in the experiment. Be- 
sides, for a < —0.5, a broad peak located at T ~ J2 in the 
temperature dependence of susceptibility emerges, which 
is a characteristic feature of one-dimensional Heiscnbcrg 
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FIG. 7: Comparison of the TMRG results for the susceptibil- 
ity in the AF-AF case with the experimental measurement in 
Ref.H 
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FIG. 9: Comparison of the TMRG results for the susceptibil- 
ity in the AF-F case with the experimental measurement in 
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FIG. 10: Comparison of the TMRG results for the magne- 
tization in the F-AF (circles) and AF-F (crosses) cases with 
the experimental measurement in Ref. ITU 



antifcrromagnctic chainiSa This suggests that the antifer- 
romagnetic component becomes significant for larger a. 
However, we failed to obtain a satisfactorily stable result 
below the peak temperature when a < — 1. 

For the AF-F case (Fig. |J, at Ji = 12.31 K and 
a = —3.5, the numerical result seems to fit the experi- 
ment data rather well, except for a little deviation in the 
left side of the peak at low temperatures. But there are 
several questionable points remained. According to the 
analysis of Mizuno et ali£, the NNN interaction through 
Cu-O-O-Cu should be antiferromagnetic, i.e. J2 > 0. 
However, the interactions of AF-F type may also be a 
candidate. In Ref. |4?J, Matsuda and coworkers proposed 
a model to explain the anomalous magnetic excitations 
in the edge-sharing CuC>2 chains of La5CagCu2404i . The 
intra- and interchain interactions are of AF-F type. (See 



also Fig. EH) 

We also calculated the field dependence of the mag- 
netization for the two cases: J\ = —138 K, a = —0.37; 
and Ji = 12.3 K, a = —3.5, to verify further if the only 
J\- J2 model is sufficient to describe the behaviors of the 
compound. We find the anomalous slow saturation of 
the magnetization cannot be reproduced either, and the 
result for the AF-F case seems worse. 

From the TMRG numerical results on the J1-J2 model 
and the comparison made with the experiment, we feel 
that only the NN and NNN interactions can not describe 
the properties of Rb2Cu2Mo30i2 satisfactorily, although 
a model based on the F-AF interactions seems to be an 
appropriate starting point. 

The above facts reveal the complexity of the interac- 
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tions underlying the edge-sharing copper oxides. On one 
hand, due to the strong electron correlations in these 
so-called Mott insulators, very limited information on 
the electronic structure can be obtained reliably. This 
makes it difficult to calculate accurately the superex- 
change interaction For instance, based on a three-band 
Hubbard Hamiltonian and cluster calculation, the NN 
and NNN interaction J±, J2 for Li2Cu02 were obtained 
as Ji = -100 K, J 2 = 62 K (a = -0.62) M The re- 
sults of the quantum chemical calculation were given as 
Jj = -142 K, J 2 = 22 K (a = -0.15)41 On the other 
hand, because the NN coupling in edge-sharing copper 
oxides is extremely small, other interactions, such as 
quantum frustrations, weak interchain correlations, and 
anisotropics can all have a chance to play an unncgligiblc 
role in determining the phase and behavior of the system. 
They are also closely related to the lattice structures and 
chemical compositions. These combined effects make any 
reasonable analysis intricate. In order to quantitatively 
recover the experimental data and various magnetic or- 
ders at low temperatures, more parameters are needed. 
This brings some "flexibility" to the theory. For in- 
stance, the broadening of the magnetic excitations found 
in Ca2Y2Cu50io requires the introduction of the antifer- 
romagnetic interchain interactions and anisotropics for 
superexchange interactions 4a In order to understand the 
helicoidal magnetic order in NaCu202, four parameters 
including frustrated longer-range exchange interactions 
are needed^! However, we wish to emphasize here that 
accurate knowledge on the behavior of model [3 espe- 
cially at low-temperature region, is indispensable in un- 
derstanding the properties of these materials. 

In conclusion, we go back to the compound 
Rb2Cu2Mo30i2 and take a closer look at the lattice 
structure of Rb2Cu2Mo30i2- Since the nearest neighbor 
Cu-Cu bond has two slightly alternating configurations 
by turns, the chain is distorted into a zigzag shape (See 
Fig. 1 in Ref. IT3) . Additional antisymmetric exchange 
interactions, such as Dzyaloshinskii-Moriya (DM4£ or 
anisotropic interactions^ together with an alternating g 
tensor, should become more important than in a straight- 
line chain. As discussed by Dzyaloshinskii and Moriya^ 
for the magnetic crystals with lower symmetries, the ef- 
fect of this antisymmetric exchange should become more 
manifest than the exchange anisotropy. The magnitude 
of this interaction is estimated as 

D~(Ag/g)J, (13) 



and the usual term of DM interaction can be expressed 
as 

X T V |S < ( 14 ) 

3 

Because of the alternating g tensor and the DM inter- 
action, an external uniform magnetic field can induce an 
effective staggered field. If D_, takes the form (— 1) J 'D, 
as it should be in the present case, the transverse com- 
ponent of the staggered field, which is perpendicular to 
the uniformly applied field, becomes dominant. For the 
S = 1/2 Heisenberg antiferromagnetic chain, a gap is 
generated by a staggered field4^ and the magnetization 
becomes gradually saturated for large fields4& The simi- 
lar mechanism may also work for the F-AF case. A com- 
bination with the frustration effect caused by the NNN 
J2 make the situation more interesting. 

Recently, it was suggested that for a < —0.38, there 
exists an incommensurate-commensurate transition at 
some critical field in the magnetization process^! The 
sharp increase up to M ~ 0.4 at B ~ 14 T and the 
following gradual saturation found in the experiment is 
argued to be connected to this phase transition. 



V. SUMMARY 

In this paper, we explored the properties of the zigzag 
spin chain with different combinations of ferro- and an- 
tiferromagnetic interactions between the NN and NNN 
sites. The existence of the gap in the F-AF case and 
the nonexistence of the gap in the AF-F case were dis- 
cussed. Thermodynamic properties of the zigzag spin 
chain in various phases were studied by using the TMRG 
method. The obtained results were used to compare with 
the experimental data from Rb2Cu2Mo30i2- We pointed 
out that for such edge-sharing copper oxide chains, be- 
sides the NN and NNN couplings, more ingredients, such 
as the interchain or DM exchange interactions, may be 
important in these materials. 
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